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Here we have studied first and second-order intertwining approach to generate isospectral partner 
potentials of position-dependent (effective) mass Schrodinger equation. The second-order intertwiner 
is constructed directly by taking it as second order linear differential operator with position depndent 
coefficients and the system of equations arising from the intertwining relationship is solved for the 
coefficients by taking an ansatz. A complete scheme for obtaining general solution is obtained which 
is valid for any arbitrary potential and mass function. The proposed technique allows us to generate 
isospectral potentials with the following spectral modifications: (i) to add new bound state(s), (ii) 
to remove bound state(s) and (iii) to leave the spectrum unaffected. To explain our findings with the 
help of an illustration, we have used point canonical transformation (PCT) to obtain the general 
solution of the position dependent mass Schrodinger equation corresponding to a potential and 
mass function. It is shown that our results are consistent with the formulation of type A A/"-fold 
supersymmetry [l4|, 0i| for the particular case AT — 1 and Af = 2 respectively. 

PACS numbers: 



I. INTRODUCTION 

There is a growing interest nowadays to design systems whose Hamiltonians have given spectral character- 
istics. In this context, the idea of designing potentials with prescribed energy spectra is worth investigating. 
Some progress in this area has been done by restricting the construction of potentials isospectral to a given 
initial one except a few energy values through the usage of Darboux transformation [l[ , factorization method 
0, supersymmetric quantum mechanics (SUSYQM) 0,|3| and other related techniques. The underlying idea 
of most of these procedures has been summarized in an algebraic scheme known as intertwining approach. In 
general, the objective of the intertwining is to construct the so-called intertwining operator C which performs 
an intertwining between an initial solvable Hamiltonian % and a new solvable one T-L with slightly modified 
spectrum such that 

CU = UC, $(x) = Ci>(x) (1) 

The ingredients to implement the intertwining are seed solutions of the initial stationary Schrodinger equation 
associated to factorization energies less than or equal to the ground state energy of H. If £ is a first order 
differential operator, the standard SUSYQM, with supercharges built of first order Darboux transformation 
operators, and the factorization method are recovered. On the other hand, if hi gher order differential operators 
are involved in the construction of C, it gives rise to higher order SUSYQM [2J]. It is possible to generate 
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families of isospectral Hamiltonians by either of the two ways: 

(i) by iteration of first order Darboux transformations. Every chain of Af first order Darboux transformation 
creates a chain of exactly solvable Hamiltonians / Hq,T-Li, / H2---T-Lm @j- Hence the intertwining operator 
between the initial Hamiltonian Ho and the final Hamiltonian %n can always be presented as a product of Af 
first order Darboux transformation operators between every two juxtaposed Hamiltonians "Ho, "Hi, ...,TIn, 

C<M = £tf£tf- 1 ...C 2 Ci, H p £ p = CpHp-!, p = 1,2, ...,Af 

(ii) by looking for the Af-th order intertwining operator directly, expressing the intertwiner as a sum of Af + 1 

terms gi(x)-^~, i = 0, 1, • • - Af, and solving the system of equations resulting from the intertwining relationship 
for the gi(xYs. 

At this point it is appropriate to mention that the quantum mechanical systems with position dependent 
(effective) mass [|| have attracted a lot of interest due to their relevance in describing the physics of many 
microstructures of current interest such as semiconductor heterostructures Q, quantum dots [8|, helium clus- 
ters and metal crystals Q etc. Recently, the intertwining operator method has been applied to Schrodinger 
equation with position dependent (effective) mass to construct first-order and chains (iterations) of first-order 
Darboux transformations and the connection between the first-order Darboux transformation and effective 
mass supersymmetry (factorization) was shown [10]. Subsequently Darboux transformations of arbitrary 
order for position dependent mass Schrodinger equation was derived and factorization of the n-th order trans- 
formation into first order transformations and existence of a reality condition for the transformed potential was 
shown 11]. In the standard supersymmetric (SUSY) approach for effective mass Hamiltonians [12j . the ladder 
operators are taken as first order differential operators similar to constant mass case, but now they depend on 
both superpotential and mass function. As a result one obtains two partner potentials with the same effective 
mass sharing identical spectra upto the zero mode of the supercharge. Second order supersymmetric approach 
(2-SUSY)was used in ref [131 ] for describing dynamics of a quantum particle with a position dependent mass. 
A compact expression for 2-SUSY isospectral pairs was derived in terms of senond order superpotential and 
the mass function. A detailed analysis has been given about zero mode equations of second order supercharges 
and possible reduction of 2-SUSY scheme to first order SUSY. Recently, a generalization of standard SUSY 
known as higher derivative SUSY or A^-fold SUSY [TiJ was given for position dependent mass Hamiltonians. 
This method keeps the basic supcralgcbra intact but differs from the standard first order SUSY in that the 
supercharges are represented as Af-ih order (Af > 1) differential operators. 

In this article an attempt is made to generate isospectral potentials of position dependent mass Schrodinger 
equation (PDMSE) by applying first and second order intertwining technique. Specifically, the second order 
intertwiner is constructed by taking it as 

Id 2 d 

where M{x) is the mass function and rj(x), ^(x) are to be determined. Substituting this in the intertwining 
relationship (JT|), we have been able to solve the apparently intricate system of equations for r](x) and 7(2;) 
by assuming an ansatz. As mentioned earlier, the closed form formulas for Af-th order intertwining operators 
and a pair of isospectral Hamiltonians with position dependent mass was already reported by Tanaka [l4j , for 
an arbitrary value of Af without recourse to any ansatz. 

The motivation for constructing second order intertwining operator directly comes from the observation 
that although an Af-th order intertwining operator can be expressed formally as a product of Af first order 
intertwining operators, it does not necessarily mean that a system constructed by an Af-th order intertwining 
operator is equivalent to one constructed by Af successive applications of each 1st order operator. In fact, it 
was shown in ref [l5j . by comparing the two approaches in the constant mass scenario, that the former is more 
general than the latter. The advantage of the direct method used here over the iterative method is that one 
can generate second-order isospectral partner potentials directly from the initial potential i.e. one need not 



3 



go through the first-order intertwining technique. We shall see that as in the case of constant mass scenario, 
it is possible to generate isospectral potentials with some spectral modification: (i) to add new bound state(s) 
(ii) to remove bound state(s) and (iii) to leave the spectrum unaffected. In this context a natural question 
is: What is the utility of finding the isospectral potentials in the position dependent mass background? To 
answer this let us note that in different areas of possible applications of low dimensional structures as already 
mentioned, there is need to have energy spectrum which is predetermined. For example, in the quantum well 
profile optimization, isospectral potentials (by deleting or creating bound states at a particular energy of the 
original potential) are generated through supersymmetric quantum mechanics. This is necessary because a 
particular effect (such as intersubband optical transitions in a quantum well) may be grossly enhanced by 
achieving the resonance conditions e.g. appropriate spacings between the most relevant states and also by 
tailoring the wave functions so that the (combinations of) matrix elements relevant for this particular effect 
are maximized [l6j . This is particularly important for studying higher order nonlinear processes. 

The organization of the paper is as follows: in section II and III we have explained the first and second-order 
intertwining techniques respectively with possible spectral modifications, with the help of a suitable example 
given in the Appendix. Also, the connection of our approach to type A Af-fo\d SUSY is shown in these sections. 
Section IV is kept for discussions and comments. 



II. 



FIRST ORDER INTERTWINING 



We consider the following two one-dimensional effective mass Schrodinger Hamiltonians (Bendaniel-Dukc 
form) [17] with the same spectrum but with different potential 



and 



Hip = Eip , H 



ntp = Eip , H = - 
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V(x) 



V(x) 
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(3) 



We connect the Hamiltonians @ and © by means of the intertwining technique. To this end, we look for an 
operator C that satisfies the relation (TTJ) . Without loss of generality let us consider the first order intertwining 
operator [Io| . as 



C = 



I 



A{x) 



y/Mixjdx 

Now using the intertwining relation (fTJ) and equating the coefficients of like order of derivatives we obtain 

2A' Ml' 2 M" 



and 



V = V 



A(V-V) 



I'M 4M 3 2Af 2 



(4) 



(5) 



A'M 1 

' M 2 



V A" 



Vm m 



(6) 



where 'prime' denotes differentiation with respect to x. 
Now using (JS|) , equation reduces to 



A" 



A'M' AM" 



Mi 



2Mi 



3AM' : 
4Mi 



2AA' + V' = 



(7) 



Integrating equation ([7]) we get 



<M 2Mi 

where fi is a constant of integration. Now we substitute 



A> AM ' ,2 rr 

-A 2 + V = li (8) 



where K=K(x) being an auxiliary function, in equation (J5J we obtain the following Riccati equation 

K' KM' K 2 , s 

Tj + V = fj, (10) 



and (fT0|) we get 



M M 2 M 



The equation (jTUJ) can be linearized by the substitution K(x) = tttj- Substituting this value of K in equations 



A(x) = -^L- (11) 



and 

-JI u "-{ai)' u ' + vu = ^ (12) 

respectively. The equation (|12l) is similar to equation @ with E = fi, fj, is sometimes called factorization 
energy and U(x) is called seed solution. It should be noted here that U(x) need not be normalizable solution 
of (|T2"j) . However, for A{x) to be well defined (without singularity), U must not have any zeroes on the real 
line. For this, we shall restrict < Eq throughout this article, E$ being the ground state energy eigenvalues of 
the equation ^ . Once we can determine the solution U of (|12l) then we shall able to construct the intertwiner 
C, the isospectral partner V and its bound state eigenvalues ip(x) with the help of following relations 

- 2U" ZU' 2 WW M" MI' 2 

MU + MU 2 + IvPU + 2M2 - 4M^ ( ' 



1 



~M \ dx U 



and 



M 



— - (InU)' 
ax 



(14) 



(15) 



The intertwiner C cannot be used to generate wave function of H at the factorization energy /it, because 
CIA = 0. We are showing below with the help of supersymmetry how ip^ can be obtained from the relation 
C^ipfi = 0, & being the adjoint of £ and is given by 
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For this we calculate & C and given by 



and 



Id 2 

n r = 

M dx 2 



M dx 2 



M' d 
Jl 2 dx' 



M' d 
APdx' 



U" M'U' \ 
~MU ~ M 2 U ) 



2U' Z 
MU 2 



U" M" 



3M r 



MU 2NP 4M 3 



(17) 



(18) 



Now from the equation ([121) we have 



V 



U" M'W 



MU M 2 U 



Substituting this value of V in (H"3l) we obtain 



V 



U" 2W M" 3M 



T'2 



MU MU 2 2M 2 4M 3 
Now using (HU) and (J2UJ) in <[TTJ) and JISJ we get 

1 d 2 M' d 



and 



M dx 2 M 2 dx 



nM 1 d 2 M' d 
M dx 2 M 2 dx 



(19) 



(20) 



(21) 



(22) 



respectively. It is clear from the equation (|22j) that the wave function of T-L at the factorization energy \jl can 
be obtained by = i.e., 



oc exp 



W M' 
~~U + 2M 



dx 



U 



(23) 



It is to be noted that if U corresponds to the bound state of H, the wave function ^^(x) defined in ([23)) is not 
normalized so that \x does not belong to the bound state spectrum of If U corresponds to the ground state 
wavefunction of % then then the potential V has no new singularity, except the singularity due to V, provided 
M is not singular and M / 0. However, if we consider U to an arbitrary state other than ground state of 
then V might contain extra singularities, which are not present in V. If U is nodeless and unbounded at 
the both end points then tjj^ix) defined in (1231) is normalizable, so that \x can be included in the bound state 
spectrum of H to generate V. In this case maximal set of bound state wavefunctions of H are given by {tp/j,, Cip}. 



A. First-order intertwining and type A 1-fold SUSY 

To show that the results obtained in the previous section are consistent with the results of type A A/"-fold 
SUSY, we are going to mention the brief results of type A A/"-fold SUSY formalism ( for details see [lH and 
references there). Type A N- fold SUSY is characterized by the type A monomial space 



9 M = {l,z,...z M - 1 ) 



(24) 



G 



preserved by Hj\f : 

d 2 „ ,/ 



n ^ = - A{z) d^~ B{z) Tz~ C{z) (25) 



where 



A(z) — a4Z 4 + (I3Z 3 + a^z 2 + a\z + do, Af > 3 
5(z) = Q(z)-^4'(z) 

C (z) = (-^-i)(.V-2) ^ (2) _ ij£_H QI[z) + R ( 26 ) 
Q(z) = b 2 z 2 + b lZ + b , Af>2 

R,ai,bi are being constants. Applying the algorithm for constructing type A Af- fold SUSY in PDM system 



l4 |. one can construct the most general form of type A A/"- fold SUSY PDM quantum systems (H,H,C) or 
equivalently {%%, Ujj-, Vtf) ■ 



P N = M{x)-^ Y[ 

where 



V ! r d „ T , . AfM'(x) Af —1 — 2k z"(x) 
— + W{x) 



i£ 

k=0 



dx v ' AM(x) 2 z'{x) 



(28) 



v+ - y +2M {im~ 2M( X ) 2 ) (29) 

W{x) = ^x—-V(xJ + 2MW (30) 



dW M (x) _ z"(x) M(x)B(z) M'{x) 
dx ~~ 2z'(x) 2z'(x) 2M{x) 



(31) 

z\x) 2 = M{x)A(z) (32) 



and the product of operators are ordered as 

N-i 

fe=0 

The solution space of the type A Hamiltonians H.^ are given by 



^ = e- w "(l,z,... )Z ")\ z=z{x) , where H# = -W u + (Af - 1) ln\z'(x)\ - - ln\M(x)\ (33) 

It is easily seen that (V — V) obtained in ([T^]) coincides with (V + — V~) in equation (|2^|) (with Af = 1) if one 
takes (comparing C with Pi) 

i ! ""W = -^' + ^M (34) 
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B. Example of first-order intertwining 

It may be emphasized that the results mentioned in section II are most general and valid for any potential 
V{x). However to illustrate the above procedure with the help of an example we shall need non-normalizable 
solutions of (|12|) corresponding to a particular mass function M{x). In Appendix A we have used point 
canonical transformation approach(PCT) to solve the equation (fT2")l . Here we are going to construct the 
isospectral partners of the following potential obtained in Appendix A (we have considered p = X = 1 for 
simplicity) 

V(x) = [(Q + & r )2 - 1] e * + e- (35 ) 



corresponding to the mass function 



M{x) = -sech 2 [\x\ (36) 



4 V 2 

The bound state solutions and eigenstates of the equation @ are given by (see Appendix A) 



_ (2n + a + S+l)n! T(n + a + S + l) \ ^ e^x ^ ( l-£ 
~ ' r(n + o-+l)r(n + * + l) J (l + e*)^ " U + e* 1 



and 



ff w = n 2 + n(g + g + 1) + ^ + ffi + ^ , n = 0,1,2,... (38) 



respectively, where 6=1 — a + + c = 1 + cr with c > | and a + b — c + | > 0. The seed solution W(x) and 
factorization energy /i are given by 

e §x / e x \ e ^~^ x ( e x 

U(x) = a- a+b+1 2F1 a,b,c, — — - +/3- a+b _ 2e+3 2F1 [a - c + 1, b - c + 1, 2 - c, 



(39) 



(a + 6+l)c c 2 

^ = -a6+^ (40) 

respectively. The asymptotic behavior of the solution U(x) given in (1391) . at both end points ±00 are given by 



resp 

m 



where 



U(x) ~ (A ia + Bx^e- ^"-^ 1 * + (A 2 a + B 2 /3)e- £ as x -> 00 (41) 



r(c)r( c - g -b) r(2- c )r( c -a-6) 

^1 = tt7 TTTTT 7 > B l = 



r(c-6)r(c-a) ' r(i-a)r(i-6) 



r(c)r(a + fe-c) r(2-c)r(a + 6-c) 

-A? = , rSo — 



and 



r( a )r(6) ' r(a-c + i)r(6-c+i) 

~ ae^-F^ e^ 1 "^ as z -> -00 (42) 
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From these asymptotic behaviors it is clear that U(x) will unbounded at x — > oo if \a + b — c\ > 1 and it is 
unbounded at x — > — oo if c < or c > 2. Therefore U(x) will nodeless at the finite part of the x axis if 
A\a + Bi(3, A 2 a + B 2 (3,a and /3 are all positive and \a + b — c\ > 1, c < or c > 2. 

Now we are going to generate isospectral potentials of the potential (13"51 with various possible spectral modi- 
fications. 

Deletion of the initial ground state : In this case the factorization energy [i is equal to the ground state 
energy Eq giving a and/or b = and U[x) becomes the ground state wavefunction i/jq(x) which is obtained 
from ([57)1 as 

U(x) = ip (x) oc (43) 

(1 + e x )^~ 

The isospectral partner of V{x) given in equation (|35)l . is obtained using equation (fT3|) . (|35)l . (f36f and (|43|) 
and is given by 

-, , c 2 - 1 _, (a + b-c)(2 + a + b-c) x a + b . A A . 

V(x) = —^-e + ± ^ >-e x + — (44) 

The above potential (|44p can also be obtained from the initial potential ([33]) by making the changes a — > 
a + l,b—>b+l,c—>c+l, this property is known as shape invariance Q. Since ipo(x) is bounded solution 

ipfiix) — is unbounded at x — > ±oo, so we have deleted the ground state energy of % to obtain V(x). 

Therefore the eigenvalues of H are given by 

E n = E n+1 = (n + l) 2 + (n + l)( Q + b)+ C(a + 6 ~ C+1) , n = 0,l,2... (45) 

Corresponding bound state wavefunctions of V(x) are obtained using equation (fT5"|) as 

K 7^ ^S^T—7^ V n' a+b ~ C+1) (-tanh^) , n = 0,l, 2... (46) 

We have plotted the potentials V(x) given in ([33]) and V(x) given in (E| for a = 5, 6 = 0, c = 3, a = 1, /3 = 
in figure 1. 




FIG. 1: Plot of the potential F(x) (solid line) given in (|35[) and its first order isospectral partners F(a;) (dashed line) 
given in (|44|l by deleting the ground state Eo = 4.5, we have considered here a = 5, 6 = 0, c = 3, a = 1, P = 0. 
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Strictly isospectral potentials : The strictly (strict in the sense that the spectrum of the initial potential 
and its isospectral potential are exactly the same) isospectral potentials can be generated with the help of 
those seed solutions which vanish at one of the ends of the x-domain. Now for f3 = and a > 0, it is seen 
from (|4*Tj) that U is unbounded at i - > oo if |a + 6 — c| > 1. But the solution ([57)1 become unbounded for 
a + b — c < — 1. So we must take a + b — c > 1. On the other hand from (|42p it is observed that U(x) — > at 
x —> — oo if c < 2 or c > but ip n (x) are not normalizable for the values of c < 2 so we must take c > 0. So 
U{x) vanishes at x — > — oo and unbounded at a;— > oo ii a + b — c>l and c > 0. In this case the spectrum of 
the isospectral potential as well as original potential are identical i.e. E n = E n , n = 0, 1, 2... Considering the 
seed solution as 

U(x) = a+b+1 2 Fi a, b, c, — — - , a + b-c>l,c>0 



(l + e*) 4 *^ > -to j 

we have calculated the explicit form of the partner potential using (fl3|) as 
1 



V(x) = ^[2c(c - 2)e- x + 2((a + b - cf - l)e x 

^c^L^ -^ 1 + c) ( 2Fi i 1 + a < 1 + M + c ' ) 2 + 4a&c 2Fi ( a ' & ' c ' 

^(o + 1)(6 + 1) (2F1 (2 + a, 2 + 6, 2 + c, ^) ) ' - (1 + c)sinhx ( 2 Fi (a + 1, 6 + 1, c + 1, ^) ) ^ 
- 8c 2 (l + c)cosh 3 x (2F1 (a,b,c, t^st)) (( a + b)cosh(^) + (l + a + b- 2c)sinh(§))}} 

(47) 

In particular for a = 3, = 5, c = 4, a = 1, /? = and using (|35[) . ([T3"l) we have obtained 

V(x) — -(23coshx + Tsinhx) 

n*) = ^ + ^ + ^# 

respectively, which are plotted in figure 2. In this case eigenfunctions and eigenvalues of the above partner 
potential V(x) are given by 

e 3x ((2 + e*)(n + 7)V { ^ (-tank (§)) + (1 + e x )(5 + 3e x )rL m {-tanh{% ))) 

UX) K (l + e^(2 + e *) se cfe(§) ( 49 ) 

£J„ = n 2 + 8n+ 10 , u = 0, Ir- 
respectively. 

Creation of a new ground state : In this case we shall consider /i < Eq. The new state can be created 
below the ground state of the initial potential with the help of those seed solutions which satisfies the following 
two conditions: (i) it should be nodeless throughout the ^-domain and (ii) it should be unbounded at both 
the end points of the domain of definition of the given potential V(x). From the asymptotic behaviors of 
the seed solution U, given in equations (|4T|) and (|42|) we have, for \a + b — c\ > 1 together with either 
c < or c > 2, the above two conditions are satisfied. But to get ip n ( x ) as physically acceptable, we shall take 
c > 2 and a + b — c > 1. In this case the spectrum of the partner potential is {[i, E n , n = 0, 1, 2...}, E n being 
the energy eigenvalues of the original potential V(x) given in (|35[) . Corresponding bound state wavefunctions 
are {^(x),^ n (x),n = 0, 1,2, ...}, where ^ and ip n are given by (f2"3"|) and (fTS]) respectively. 
For a + b — c>l,c>2 and the seed solution U given in (|3"9")) , the general expression of the isospectral potential 
becomes too involved so instead of giving the explicit expression of the partner potential we have plotted in 
figure 3 the original potential V(x) given in (|35[) and its partner potential V(x) (which is obtained using (fT3|) ) 
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-4 -2 2 4 

FIG. 2: Plot of the potential V(x) (solid line) and its first order isospectral partner (dashed line) given in (|48l) . 

considering the particular values a — 2.8, b — 20, c = 4.4 and a = /3 = 1. In this case the energy eigenvalues 
of V(x) are given by 

E n = {-13.32, E n , n = 0, 1,2...} = {-13.32, n 2 + 22.8n + 42.68, n = 0,l, 2...} (50) 
Corresponding eigenfunctions can be obtained using the formulae (|23[) and (|15D . 




-4 \ -2' 2 4 6 



FIG. 3: Plot of the potential V(x) (solid line) given in (|35[l and its first order isospectral partner V (dashed line) by 
inserting the state [i, = —13.32. We have considered here a = 2.8, b = 20, c = 4.4, a = /3 = 1. 



III. SECOND ORDER INTERTWINING 



Now we assume the existence of a second order intertwining operator 



Id 2 d 



where rj(x),^/(x) are to be determined. Substitution of this intertwiner in equation ([T]) and comparison 
coefficients of like order derivatives leads to a set of following equations 

, M' 3M' 2 2M" 
V = V + 2rf + —r) 



M ' M 3 M 2 
2V 2 7 ' rf' rfM' M"tj 2M' 2 r] 6M' 3 QM'M" M'' 





tj'M' 


M"n 
f ' 


2M' 2 r) 


M 


M 2 


M 2 


M 3 


7" 


M'-y ' 






M 


M 2 







M M M M 2 M 2 M 3 M 5 M 4 M 3 



M 



Now using (I52j) the equations (153)) and (|54[) reads 

, M'rj 2 3r]M' 2 2M"r) 2 7 ' M'r( 

m M 3 " + ~W 2 W + ~W 

2M' 2 j) rf' M"r) 2V 6M' 3 SM'M" M'" 

+ M 3 M ~ M 2 ~ ~M~ ~ M 5 + M 4 ll 3 ~ 

and 

/ , M'n 3M' 2 2M"\ M'y' 7" . V" 
7 27,' + —L-— - + —-)+ —J- - J- - r,V - — = 



M M 3 M 2 J M 2 M 1 M 

respectively. Equation (|55|) can be integrated to obtain 

_ Mr] 2 M't] rf M' 2 M" 

7 ~ ~Y~ + ~2M ~ ¥ ~ M 3 " ~ 2M 2 + 1 

where C\ is an arbitrary constant. Using (|57[) in (|56l) we obtain 

rf" i j r , 2 rjri'M' , 2 n , T , , 5r/M' 2 3r/M" , r; 3 M' 7y 2 M' 2 #M' 2t?M /: 



2M +M? ^ 2M 2?/ ' 2??V ' M 3 M 2 2 2M 2 M M 4 

rfM" 5r]M'M" rf'M 1 .. . r\M"' . C^M' WM' 2 18M' 4 

+ ^M- + ^M 3 --^P- T]f] ~ 2NP +2Clfl + ^T + ^T 3 M6- 

49M' 2 M" 3CiAf' 2 2VM" AM" 2 2C X M" V'M' 9M'M"' M"" _ 

+ 2M 1, M 3 M 2 M 4 "" + M 2 M 2 2M 4 + 2M 3 ~ 

Multiplying by -TT- ) , above equation (|58[) can be integrated to obtain 

Mr) V -\- U-iMr) -\ 5 3 

4 M 3 M M M 3 2M 

!\/TW ^2^2 M y/ M /y M /» r/ 4M'M"r) 3M' 4 5M' 2 M" 



7777 tT 




2 4 


- - »/V • 


M'rf 


5M' 3 r] 


+ 2 4 


m 4 


M" 2 


M'M'" 


4M 4 + 


2M 4 



+ 



M AM 2 2M 2 2M 2 2M 2 M 3 M 6 2M 5 

C 2 = 



12 



where C 2 is the constant of integration. For a given potential V(x), the new potential V(x) and j(x) can be 
obtained from (1521) and (|57|) if the solution 7y(a;) of (|59|) is known. To obtain t](x) we take the Ansatz 

2 / M'\ M' 2M' 2 M" , , . 

^'^V + 2(, + ^jr + -, + ^- — + ^ (60) 

where £ is a constant to be determined and r is a function of x. Using above ansatz in equation (|59[) we obtain 
the following equation 



t' , t 2 M't ta { ^ ^ 2 , 2M ' ^ r2 M ' T T r , ^ £ 



M 3 \M M M 2 2 \ 4 



(61) 



Since equation (|61[) is valid for arbitrary 77, the coefficients of each power of 77 must vanish, which give £ 2 = 4C*2 
and 

M M M 2 2 v ; 

Now defining /j = C\ — § , the above equation can be written as 



T T 



2 



M't jr £ 



= V — fi , it = d - ^ = Ci T VC2 (63) 



M M M 2 ^ ' A 2 

The equation (|63p is a Riccati equation which can be linearized by defining t = Making this change in 
equation (|63|) we obtain 

- ^W" - (^) U' + VU = iU (64) 

Depending on whether C2 is zero or not, £ vanishes or takes two different values ±yC2- If C2 = 0, we need 
to solve one equation of the form (|f)3"j) and then the equation (jSU|) for r?(a;). If C2 ^ 0, there will be two 
different equations of type (|63| for two factorization energies /ii,2 = C\ =F n/C^- Once we solve them, it is 
possible to construct algebraically a common solution rj(x) of the corresponding pair of equations (|60|). There 
is an obvious difference between the real case with C2 > and the complex case C2 < 0; thus there follows a 
natural scheme of classification for the solutions rj(x) based on the sign of G%. In our present article we shall 
not discuss the case C 2 = 0. 

(i) Real Case (C 2 > 0) 

Here we have /ti^ € R, Hi ^= [i 2 - Let the corresponding solutions of the Riccati equation (|63p be denoted by 
Ti :2 (x). Now the associated pair of equations (|60l) become 

9 / M'\ M' 2M' 2 M" , . 

v > = Mr, 2 + 2l V+1 ^)T 1 + — 11+ — - — + ^ 2 - lll (65) 



and 



, M'\ M' 2M' 2 M" 
MV 2 + 2(V + JI , U + _„ + _-_ + Ml - M2 (66) 
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respectively. Subtracting (|55|) from and using we obtain ry(x) as 

^ W ~ri-r a M 2_ MW(U 1: U 2 ) ( ' 

where U\,Ui are the seed solutions of the equation (|64[) corresponding to the factorization energy [i\ and /i2 
respectively and W{U\,U2) — UxU'2 — U[U2, is the Wronskian of U\ and U2- 

Now it is clear from ([521 and (|67l) that mass function M(x) is nonsingular and does not vanish at the finite 
part of the x-domain, so that the new potential V(x) has no extra singularities (i.e. the number of singularities 
in V and V remains the same) if W(U\,U2) is nodeless there. The spectrum of % depends on whether or not 
its two eigenfunctions ij)^ 2 which belongs as well to the kernel of C) can be normalized [21| . namely 

tftfto = and Hrf^ = fij^ H , j = 1,2 
where & is the adjoint of C and is given by 

M dx 2 \' M 2 J dx \ M 3 M 2 1 7 

For j = 1 the explicit expression of the two equation mentioned in (|68p are 

1 d 2 ^ ( 2M'\d^ , (2M> 2 M" \ 

■ V +7 h/Vi = ( 69 ) 



M dx 2 V /da; V M 3 M 2 

and 



1 / 1 

M M1 ~\M 



C-lXf ) ^+(^-^1)^=0 (70) 



respectively. Adding (p9")) from ((70)1 we obtain 

Af \ dib u , ( - 2M' 2 . M"\ , 

Substituting the values of V" and 7 from (f52"j) and (|57| with 2Ci = /ii + /12, in the above equation (|TTj) . we get 

dz ( '° 5VW 2(„ + ^) (72) 

Now using our ansatz (f^Uj) in (|T2"|) and then integrating we obtain 

^ K — is — K "^m^y (73) 

Above procedure can be applied to obtain as 

7 tjM + MUi 

ib u , oc sL. oc — ±_ (74) 

If both 2 are normalizable then we get the maximal set of eigenfunctions of 1-L as{VVii j VV2 j VVi oc Cif> n }- 
Among the several spectral modifications which can be achieved through the real second order SUSYQM for 
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PDMSE, some cases are worth to be mentioned. 



Deletion of first two energy levels : For — Eq and ^2 = Ei the two solutions of equation 
are the normalizable solutions of equation ^ i.e, U\ = ipo(%) and U2 = ipi(x) respectively. It 
turns out that the Wronskian is nodeless but two solutions VVi anc ^ V^a are non-normalizable. Thus 
Spifi) = Sp(H) — {E a , Ei} — {E 2 , E3, £4, ...}, i.e., the two levels E and E\ are deleted to generate V. 

Isospectral transformations : If we take [i\ < [12 < Eq and choose U\ and IA2 such way that either 
Ui,2{x{) = or 1Ai t 2{x r ) = 0, xi and x r being the end points of the domain of definition of V(x), then 
the Wronskian W(IA\,IA2) vanishes at xi or x r . Hence an< ^ VV2 become non-normalizable so that 
Sp(H) = S P (H). 

Creation of two new levels below the ground state : For n% < A*i < Eq and choosing U\ and 
U2 in such way that U2 has exactly one node and U\ is nodeless then the Wronskian W(JA\,Uz) becomes 
nodeless, also two wavefunctions ip^ and ip^ 2 are normalizable. Therefore the spectrum of H becomes 
S'p('H) = Sp{'H) ^2} = {fJ-i-, (J-2, E n , n = 0, 1, 2...} i.e. two new levels have been inserted to the spectrum 

of V(x) to obtain V(x). 



(ii) Complex case (C2 < 0) 

For C2 < the two factorization energies fi\ and /12 become complex. In order to construct real V we shall 
choose /ii and fi2 as complex conjugate to each other i.e, [A\ = fx € C and [12 = p.. For the same reason we 
shall take ti(x) = t(x) and T2(x) = f(x). Hence the real solution i](x) of (|60|) generated from the complex 
t(x) of (1551) becomes 

_ /i - /2 M< _ Jm( M ) M' W^g) 

— 7 = ~ 772 ~ 77777 772^ 777777777 V ^ 1 



M 2 Im{r) M 2 MW(U,U) 

__________ ml K 121 r>rM"»n oc 

M(ju-m) 



Defining = Yif^-X ' becomes 



»»(*) = "777 - 772 (76) 



Mw M 2 

For the factorization energies /i and /2 the equation (|64|) becomes 

— l -U" - (—} U' + VU=u,U and - — U" - ( — ) W + VU = fiU 
M \M I M \ M 



Multiplying first equation by U and second equation by U and then subtracting we obtain 

= l«| 2 (77) 



W'{U,U) M'W(U,U) _ , r i 



M(n-fl) M 2 ^-/2) 
Using above relation (|TT[) we have 



which implies that w(x) is a non-decreasing function. So it is sufficient to choose 

lim U = or lim U = (79) 

for the Wronskian W to be nodeless. It is to be noted here that in this case we can only construct potentials 
which are strictly isospectral with the initial potential. 
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A. Second-order intertwining and type A 2-fold SUSY 

The second-order intertwiner C in equation (IBTj) coincides with P2 given in equation (|28|) if one takes 

2W(x) M'(x) 

*» (<B) = MM " MW (80) 
It is now easy to verify that for this rj(x), (V + — V~) given in ((29)) (with Af = 2) agree with (V — V) given in 

Now it is to be shown that the Hamiltonian "H given in equation ^ admits two eigenfunctions Ui2{x) 
corresponding to two factorization energies fj,% t 2 respectively i.e., 

H Ui(x) = m Ui(x), i = l,2 (81) 

will belong to type A 2-fold SUSY in PDM background (in constant mass scenario this was already proved in 
ref.[l8j]). For this we define 

= ^44, W 2 -(z) = W(z) = -In U x {x) (82) 

For this W(z), it is evident that the gauged Hamiltonian l-L^ defined by 

H2 = e w He~ w (83) 

must be diagonal in the basis V2 — (1, z) because of the assumption (|8Tj) and the choice ([82]) . From equation 
(f3T|) . its immediate consequence (for M = 2) is 



z"(x) 2f'(a;)M'(a;) 2z'{x) 2 dW{z) 
M{x) M(x) 2 M(x) dz 



B y z > = 1Z ( 84 ) 



Using equations (|82|) and (1811) in the above equation (|84)) it can be shown that 

B(z) = (pi - /i 2 ) z{x) (85) 

For this value of B(z) it is also easy to verify that the expression of r)(x) in equation (|5H1) and (pT)) are same. 
Now it is evident that the gauged Hamiltonian B.^ preserves the vector space V2 = (1, z). Hence it is possible 
to get type A 2-fold SUSY system (%,%,£) following the prescription given in ref. [14], with the choice of 
z(x), W(z) and H2 given by (|82]) and ([25} respectively. 



B. Example of second-order intertwining for real factorization energies 

It may be emphasized that the results mentioned in section III are most general and valid for any potential 
V(x). However to illustrate the above procedure with the help of an example we shall need non-normalizable 
solutions of (l64l) (which is similar to equation ([T2|) but with two factorization energies) corresponding to a 
particular mass function M(x). To illustrate the second order intertwining with an example we have considered 
the potential (|35|) as an initial potential. Corresponding seed solution for the factorization energy // = [i\ which 
is obtained in Appendix A, is 

. (a + b+l)c c 2 ,„„. 
u 1 = -ab+± >--— (86) 
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Ui( x ) = a ~ " , a+b+i 2^1 ( a,b,c,— — - + /3- „ +i ,_2c+3 2-Fi a-e+l,6-c + l,2-c, 



(1 + e 1 ) 2 V 1 t e / (1 + e^j 2 \ J- 



(87) 

We notice that the potential (|35|) and corresponding Hamiltonian are invariant under the transformation 
a — s- a + v and b ^ b — v. z/ e K. — {0}. But the solution (IMl) of the corresponding Schrodinger equation 
changes to 

B^ x ( C 

£/ 2 (a:) = a a+b+1 2 F\ a + u,b— v,c. 



l + e x 

('-§) 

+b-2c+3 



/3 a+b _ 2e+3 2-F1 I a + v — c + 1,6 — v — c+1,2 — c, 



and the corresponding factorization energy is given by 

(a + b + l)c c 2 9 
^ 2 = -a6+^ — + u{a-b) + v 2 (89) 

Thus the general solutions of the equation (|64|) for the two factorization energies fix and /12 , are given by (|87|) 
and (|88p respectivelv. The asymptotic behaviors of the seed solution U2 remains same as 1A\, which are given 
in (glj) and (H2l . 

Deletion of first two energy levels : Let us take /ii = Eq and fj,2 = E\, U\ = ^po(x) and U2 = ^i{x) which 
are given in ([37)l . The Wronskian W{U\ ) U2) is given by 

wiUM « (1+ e( e C ^ b+3 (90) 

which is nodeless and bounded in (—00, 00) as c > —\ and a + b — c + i >0 (these conditions are mentioned 
at the end of the Appendix A). The second-order SUSY partner of V(x) is obtained using equation (|52")l and 
is given by 

V(x) = J [(c 2 - 2c(a + b + 2) + {a + b+ l)(a + 6 + 3)) e x + c(c + 2) e - a + 4(a + 6 + 1)] (91) 

Clearly the eigenfunctions oc ;y(^^ 2 ) and V> M2 oc w ^^ 2 ) of "H associated to fi\ — Eq and /i2 = E\ are 
not normalizable since 

lim ip^ ia (x) = 00 



X— ¥ — 00,00 



Thus Sp(H) = - {£ ,£i} - {E 2 ,E 3 ...}. 

In particular taking a = 5, = 0, c = 3 we have plotted the potential U(x) and its second-order SUSY partner 
V{x) given in ([3"5]) and (f^Tj) respectively, in figure 4. 

Strictly isospectral potentials : The strictly isospectral partner potentials can be constructed by creating 
two new energy levels in the limit when each seed vanishes at one of the ends of the x-domain. Now from 
the asymptotic behaviors of the seed solutions, we note that both the seed solutions vanish at x — > —00 for 
/? = 0, a>0ifa + 6 — c > 1 and c > 0. Considering j3 = 0, a = 1 in ([57)1 and ([8"8"]) we take two seed solution as 

UX[X)= (l + e")W 2J?1 ( a ' 6 ' C 'TT^) (92) 



and 



(1+ 5yi±i±r 2Fl ( a + "' & -^ c 'r+-^ 



^(a) = 7T~ — ^ a+b+i ■jF l [a + v,b-v,c, - ■ - ) (93) 
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FIG. 4: Plot of the original potential (solid line) for a — 5,b = 0, c = 3 and its first-order SUSY partner (dashed line) 
by deleting the ground state Eo = 4.5 and second-order SUSY partner (dotted line) by deleting two successive states 
E = 4.5, E 1 = 10.5. 

Since Ui^{x) — > at x — > — oo, from the expressions ([751) and ([74")) we can conclude that 

lim t/j,^ (x) = oo 

x— > — oo 

which implies that /J12 does not belongs to Sp(H) i.e. V"(x) is strictly isospectral to V(x). Here the general 
expression of the partner potential is too involved so instead of giving the explicit expression we have considered 
particular values a = 3, b = 5, c = 4, v = 1, a = l,j3 = 0. Corresponding expression of the partner potential 
and its energy spectrum are 

3 

F(x) = 1 + -{Qcoshx- Isinhx), E n = E n = n 2 + 8n + 10, n = 0, Ir- 
respectively. In figure 5, we have plotted the initial potential, its first and second-order strictly isospectral 
partner potentials for the parameter values a = 3, b = 5, c = 4, a = l,/3 = 0, v = 1. 




-4 -2 2 4 



FIG. 5: Plot of the original potential (solid line) and its first-order (dashed line) and second-order (dotted line) SUSY 
partner by making the isospectral transformation for a = 3, b = 5, c = 4, a = 1, P = 0, v = 1. 
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Creation of two new levels below the ground state : Two energy levels can be created taking 
M2 < Ml < and using those seed solutions U\ and U2 for which the Wronskian become nodeless. In this case 
the expressions of the Wronskian contains several Hypergeometric function, so it is very difficult to mention 
the range of a,b,c and v for which it is nodeless. In particular for a — 2.8, b — 20, c = 4.4, a = 1,/? = 1 we 
have the Wronskian is found to be nodeless. For the same values of a, 6, c we have plotted the potential and 
its second order partner in figure 6. The second-order isospectral partner is obtained using equation (1521) . 
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FIG. 6: Plot of the original potential (solid line) for a — 2.8, b — 20, c = 4.4 and its first-order SUSY partner (dashed 
line) by creating a new level /ii = —13.32 and second-order SUSY partner (dotted line) by creating two new levels 
fii = -13.32,^2 = -85.32 
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C. Example of second order intertwining for complex factorization energies 

As mentioned earlier, in this case we can only construct the strictly isospectral partner potentials. The 
complex factorization energy fix and /i2 given by equation (1861) and (|89[) . can be made conjugate to each other 
in several ways. One of the way is by making following restrictions on a, b, c, v : c € R, Im(a) — —Im(b), v — 
Re(b) — Re(a). But in order to keep the initial potential real we have to made two more restrictions e.g. 
Re(a) + Re(b) — c > 1 and c > 2. In particular taking a = 6.1 — 5 i,b ^ 8 + 5 i,c — 4.1, v = 1.9 we have two 
factorization energy /ii(= /i) = —51.25 + 9.5 i and /ia( = P) = —51.25 — 9.5 i. For these values of a, b, c, v and 
a = 1, /5 = the seed solution W becomes 

U (*) = (1 + e , )7 . 55 ^ (6.1 - 5 i, 8 + 5 *, 41, J (94) 

Clearly U(—oo) = and — > oo as x — ^ oo so this seed U and its conjugate U can be used to obtain the 
second-order SUSY partner potential V(x) with the help of equations (|52|) and (|75|) . In figure 7 we have plotted 
the initial potential V(x) given in (|35p and its isospectral partner V(x) for the parameter values mentioned 
earlier. 




-4 -2 2 4 



FIG. 7: Plot of the original potential (solid line) for a = 6.1 — 5 i,b = 8 + 5 i,c — 4.1, u = 1.9, /ii = — 51.25 + 9.5 i, /i 2 = 
—51.25 — 9.5 i. and its second-order isospectral partner (dashed line). 



IV. SUMMARY AND OUTLOOK 

In this article we have discussed the possibilities for designing quantum spectra of position dependent 
mass Hamiltonians offered by the intertwining technique. For doing this, we start with the non-normalizable 
solution of position dependent mass Schrddinger equation with the initial potential (obtained by using the 
point canonical transformation approach). To generate spectral modifications by first order intertwining, we 
have used solutions to the position dependent mass Schrddinger equation corresponding to factorization energy 
(not belonging to the physical spectrum of the initial problem) less than or equal the ground state energy 
in order to avoid singularity in the isospectral partner potential provided the mass function is not singular 
and is not equal to zero in the real line. Thus it is possible to generate isospectral partner potentials (a) 
with the ground state of the original potential deleted (b) with a new state created below the ground state 
of the original potential (c) with the spectrum of the original potential unaffected. In ref |10| . the first order 
intertwining technique was illustrated by considering the free particle case. 
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In the case of second order intertwining, instead of using the iterative method used in [10j |. the second 
order intertwiner is constructed directly by taking it as second order linear differential operator with unknown 
coefficients which are functions of x. The main advantage of this construction is that one can generate second- 
order isospectral partner potentials directly from the initial potential without generating first-order partner 
potentials. The apparently intricate system of equations arising from the intertwining relationship is solved for 
the coefficients by taking an ansatz. In this case the spectral modifications are done by taking appropriately 
chosen factorization energies which may be real or complex. For real unequal factorization energies, it is 
possible to generate potentials (a) with deletion of first two energy levels (b) with two new levels embedded 
below the ground state of the original potential (c) with identical spectrum as of the original potential. For 
complex factorization energies, it is shown how to obtain strictly isospectral potentials. It must be mentioned 
here that in all the above cases the conditions for having spectral modifications remain the same as in the 
case of constant mass scenario [2l| provided the mass function M is nonsingular and is not equal to zero in 
the finite part of the real line. 

In this article, the equivalence of our formalism to type A AT- fold (J\f =1,2) SUSY in PDM background 
is shown. Also, it is shown that an arbitrary one body quantum PDM Hamiltonian which admits two eigen- 
functions in closed form belongs to type A 2-fold SUSY as was previously done in constant mass scenario 

[Ill- 
Some of the interesting issues to be investigated in future are 

(i) to obtain spectral changes that appear above the ground state energy of the initial potential. Specifically, 
how to create/delete a pair of levels between any two neighboring initial ones, how to move an arbitrary level 
or delete an arbitrary level. Specially interesting will be the possibility of embedding a single level at any 
arbitrary position. 

(ii) to obtain spectral modifications when the two factorization energies are equal. 



Appendix A: Construction of exactly solvable effective potential via PCT 

In order to find the general (unbounded) solution of the equation (fT2"|) we shall use PCT method in PDM 
background [13] to solve this equations. Let us find the solution of equation (|12p of the form 

U{x)=f{x) F(a,b,c,g) (Al) 

where f(x), g(x) are two function of x to be determined and F(a, 6, c, g) is the Hypergeometric function which 
satisfies second order differential equation of the type 



d 2 F , x dF . . , N c-(a + b+l)g , s ab 

- y +g J - + %F = 0, with Q(g) = ^- and R (g = -— r 

dg 2 dg 9(1-9) 9(1- 9) 



Substituting equation (|A1[) in (|12|) we obtain 

d 2 F ( g" 2f M' \ dF ( f" . M M'f 



dg 2 \g' 2 fg< Mg> ) dg \fg' 2 v ' ' g' 2 Mfg>- 
Comparing equation (|A2|) and (|A3j) we get 

Q{ 9 {x)) = £ + 

After simplification of the above equation (|A4I) we obtain 



+ (p- V )^-T7TZo) F = G (A3) 



(A4) 



/(,■) x J 9- exp U ] Q(t)dt\ (A5) 
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and 



9 



12 



ff 

M 



P V 2Mg' AM\g'J + M\ 2 dq 4/ 2A/~ 4.U ; 



ldQ Q 2 \ M" 3M' : 



2 



(A6) 



respectively. Now in PCT approach there are many options for choosing M(x) [22j |. for example M(ir) = 
Xg l2 (x) , M = Xg'(x), M = , A being a constant. Here we choose M(x) = Xg'(x). For this choice of the 



9'{x) 

mass function and using the values of Q(g), R(g) given in equation (|A2j) . equation (|A6|) reduces to 



H-V 



A 



a 6 



(c- (a + b + l)g) a + 6+1 c- (a + b+l)g c-(a + b + l)g 



g(l-g) 4g 2 (l-g) 2 2g(l-g) 2g*{l - g) 



5(1 -9? 



(A7) 



Now in order to generate a constant term on the right hand side of the above equation which will correspond 
to p on the left-hand side, we set xg(i-g) = P> wnere P is a positive constant. This gives 



9(x) 



and Mix) = sech 2 

l + eP Xx 4 



oo < x < oo 



(A8) 



For these values of g{x) and M(x) we obtain from equation (IA7|) new potential V(ie) and factorization energy 
/i as 



^ = [(a + b~c) 2 ~l]P cPXx + cp(c-2) c _ pXx ^ 



OO < X < oo 



(A9) 



and 



/! = —abp - 



(a + 6 + l)cp c 2 p 



(A10) 



Also from (|A5j) we get 



/(*) 



e 2 3 



(1 + eJ> Xx )^— 



Hence the solution of the equation (jl2j) at the factorization energy /i is given by 

p p\x 



U(x) 



e 2 



.... 2-F] a, 6, c, ; 

Another linearly independent solution of (|T2"j) at the same factorization energy can be written as [l9| 

pAan 



(All) 



a+t ,„ 2c+3 2F1 ( a-c+l,6-c+l,2-c, 



(1 + eP Aa; ) 5- 



1 + eP Aa 



The linear combination of above two solutions can be taken as the most general non-normalizable solution of 
the equation (|12p at the factorization energy p,, and is 



U{x) = a 



e 2 3 



(1 + eP Aa: ) £ 



— [a,b,c 



a p\x 



1 + P+ Al 



(1 + p+ A:r ) £ 



a pXx 



— 2F1 la — c+1,6 — c+1,2 — c, — 



a pA:r 



(A12) 
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where a and /3 are two arbitrary constants. Consequently the bound state solutions of the equation @ for the 
potential (|A9|) . are obtained from equation (|A12I) by putting a = 1, /3 = and a = — n, b = I — a + a + 5, c — 
l+a (see 15.4.6 of ref. H) 



V r(n + (T + l)r(n + £ + l) / (i + epAx)**** 3 " Vl + e pA " 
and the energy eigenvalues are given by 

E„ = n 2 p + np(a + 5 + 1) + (g + ^ + ^ , n = 0,l,2,... (A14) 

It should be mentioned here that for ip n (x) to be a physically acceptable solution it should satisfy the following 
two conditions: 

(i) It should be square integrable over domain of definition D of M(x) and tp(x) i.e., 



\tp n (x)\ 2 dx < OO 

(ii) The Hcrmiticity of the Hamiltonian ([2]) in the Hilbert space spanned by the eigenfunctions of the potential 
V(x) is ensured by the following extra condition [23] 



at the end points of the interval where V(x) and tp n (x) are defined. This condition imposes an additional 
restriction whenever the mass function M(x) vanishes at any one or both the end points of T>. In order 
to satisfy this two conditions we have to impose a restriction a > — | and S > — | or equivalently c > 
\ and a + 6- c+ | > 0. 
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